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The non-local properties of the noisy three-qubit Greenberger-Horne-Zeilinger (GHZ) states, p v = 
v\GHZ) (GHZ\ + (1 — v) 11/8 parameterized by the visibility < v < 1 are investigated. Based on the 
violation of the 2 x 2 x 2-setting Mermin inequality, p v is non-local for the parameter range 1/2 < v < 1. 
It has been posed whether additional settings would allow to lower the threshold visibility. Here we report on 
Bell inequalities giving a threshold value smaller than v — 1/2. This rules out the possibility of a local hidden 
variable model in the limit of v — 1/2. In particular, the lowest threshold visibility we found is v — 0.496057, 
attainable with 5x5x5 settings, whereas the most economical one in number of settings corresponds to 
3x3x4 settings. The method which enabled us to obtain these results, and in particular the about 10000 tight 
Bell inequalities giving v < 1/2 are also discussed in detail. 



I. INTRODUCTION 

Quantum theory allows correlations between remote sys- 
tems, which are fundamentally different from classical cor- 
relations. Quantum entanglement is in the heart of this 
phenomenon |Q]]. In particular, separated observers may 
carry out local measurements on entangled quantum states 
on such a way that the correlations they generate are out- 
side the set of common cause correlations. This is known 
as Bell's theorem yj]. 

Moreover, such quantum correlations may find applica- 
tion in novel device-independent information tasks. They 
enable perfect security IB Si or randomness generation J5[] 
without the need to trust the internal working of the devices. 

It is desirable to find states which tolerate the largest 
noise in order these protocols be useful. In this regard, let 
us note a link between the strength of violation of certain 
multipartite Bell inequalities and the security of quantum 
communication protocols @|. 

For the simplest case of a Werner state fj], which is the 
maximally entangled singlet state mixed with white noise, 
the CHSH Jst] seems to be the best inequality for a modest 
number of settings (< 10) HHU], giving the critical vis- 
ibility v = 1/ v2- Recently, this value has been slightly 
overcome with 465 number of settings floll . 

When we turn to a maximally entangled state shared by 
more than two parties, the threshold visibility drops down 
rapidly JTlll (exponentially in the number of qubits). Even 
for three parties by sharing the noisy 3-qubit GHZ state the 
critical visibility becomes v cr u — 1/2 using the Mermin 
inequality (note that the separability limit corresponds to 
v sep = 1/5 HH]). To the best of our knowledge, there 
has been no known Bell inequality giving a lower thresh- 
old. Let us note that in recent studies lfl3l [Till using a so- 
phisticated numerical method the lowest value found up to 
settings 5x5x5, matched v = 1/2 associated with the 
Mermin inequality. In the present paper, we provide a sys- 
tematic study to explore Bell inequalities with violations 



beyond the Mermin value. Our findings may also bear rel- 
evance to the problem of simulating GHZ-correlations with 
classical communication lfl5i[l6ll . 

We start by introducing the Mermin inequality (Sec. HJ). 
Then, in Sec. [Ill] a heuristic method is given to generate 
tight Bell inequalities with potentially large violations. In 
Sec. lIVl we present the results by listing the Bell inequali- 
ties according to different aspects. In Sec. [V) a symmetric 
one is highlighted, with small coefficients, and a closed so- 
lution is given for the maximum quantum violation yield- 
ing a critical visibility lower than 1/2. The paper ends in 
Sec.jVjby summarizing the results and posing open ques- 
tions. 



II. MERMIN INEQUALITY AS A CASE STUDY 

Let us represent the 2x2x2 setting Mermin polyno- 
mial ifllll in terms of three-party correlators, 

M =A B C - AqBxCx - AxBqCx - A&Co, (1) 

where the above polynomial is to be understood as a sum 
of expectation values; for example, AiBiCq denotes the 
average value of the product of the outcomes of Alice, Bob 
and Cecil when both Alice and Bob performs the second 
measurement and Cecil performs the first one. The local 
bound of this inequality is C = 2. 

For a given number of measurement settings and out- 
comes the correlations which can be described by a local 
realistic model define a polytope (the so-called Bell poly- 
tope), which is a convex set with a finite number of extreme 
points. Bell inequalities define the limits on these correla- 
tions. The Mermin inequality is a facet of the 2x2x2 
local 3-party-correlation polytope fl7ll (involving only full- 
tripartite correlation terms) and also for the 2x2x2 Bell 
polytope lfl8ll (involving single party marginals and 2-party 
correlators as well). Note that in order to get the lowest crit- 
ical visibility for a given number of settings and outcomes, 
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it is always sufficient to consider the facets of the local poly- 
tope. We refer to inequalities associated with these facets as 
tight Bell inequalities (hence Mermin inequality is a tight 
one). Furthermore, it suffices to focus on inequivalent Bell 
inequalities, that is, on those which are not equivalent un- 
der relabelling of measurement settings, outcomes and ex- 
change of parties. We would also like to mention that the 
inequivalent Bell inequalities to be presented turn out to be 
tight both in the full-tripartite-correlation space and in the 
full probability space. 

Let us now take the noisy 3-qubit GHZ state, 

Pv =v\GHZ)(GHZ\ + (l-v)l, (2) 

o 

where '\GHZ) = (|000) + |111>) /V2 is the 3-qubit GHZ 
state |19j|, and v is the visibility parameter. With equato- 
rial von Neumann measurements, corresponding to trace- 
less observables, the maximum quantum value saturates the 
algebraic limit (Q = 4). This is the no-signalling limit 
as well, hence Mermin inequality exhibits pseudotelepathy 
i2(ill . Note, that the inequalities we present in the following 
will not have this property, that is, the no-signalling bound 
will never be saturated. We denote the threshold visibil- 
ity for which quantum correlations turn to classical by v cr u 
(namely, for v > v cr i t in Eq. (f2| there do not exist local 
realistic correlations reproducing the quantum ones). From 
the form of Eq. (0 it follows the ratio v cr i t = C/Q, where 
C is the local bound and Q is the quantum bound achiev- 
able with traceless observables and yields v cr i t = 1/2 for 
the 3-qubit GHZ state based on the Mermin inequality. In 
the following we will exhibit several inequalities which go 
beyond this bound. 



III. THE METHOD 

To find Bell inequalities with low classical per quantum 
values (hence low critical visibilities), we applied a two 
step procedure. In this paragraph we describe the two steps 
briefly, and in the next ones a more detailed description is 
given. In the first step we used a simplex downhill algo- 
rithm [21] to minimize the classical per quantum value for 
Bell inequalities of very special form, characterized by just 
one parameter for each measurement setting. The initial 
parameters were chosen randomly. This procedure leads 
to non-tight Bell inequalities (that is, which do not define 
facets of the local full-tripartite correlation polytope) with 
not exceptionally good critical visibilities, usually between 
0.5 and 0.6. In no case we got a desired value of less then 
0.5 in this way. However, the best critical visibility can 
always be achieved with a tight Bell inequality, therefore, 
in the second step we determined the facet of the poly- 
tope crossed by the ray pointing towards the direction in 
the space of correlations defined by the Bell coefficients we 



have got in the first step. The tight inequality correspond- 
ing to this facet always gives a critical visibility value better 
than the original inequality does. Here we worked in the 
restricted space of three-party correlations, but all inequal- 
ities we have got this way are tight inequalities in the full 
space as well, which includes two-party correlation terms 
and single-party marginal terms. 

The coefficients of the special Bell inequalities, which 
multiply the corresponding three-party correlators AiBjCk 
we considered in the first step are given as: 

M ijk = cos^f + <pf + tp°). (3) 

When optimizing the critical visibility we have taken 
Sl=i £j=i YlT=i M i 3 k as tne quantum value of the Bell 
expression, that is the coordinates of the point in the space 
of correlations for which we calculated the quantum value 
are just the Bell coefficients themselves. This point does be- 
long to the set of quantum correlations, as it can be achieved 
by von Neumann projective measurements performed on 
components of a shared 3-qubit GHZ state. If we choose 
the ith, jth and fcth measurement operator A;, Bj and Ck 
of Alice, Bob and Cecil, respectively as 



A, 


= cos (pfa x - 


h sin ipf a y 




= cos iff a x - 


\- sin tpf o y 


c k 


= COS (fk&x - 


hsin <f^a y , 



where a x and a y are Pauli operators, then the joint correla- 
tion is given by the right hand side of Eq. (f3]), indeed, that 
is: 

AiBjCk = (Ai ® Bj ® C k ) = cos(<p? + <pf + tp%). (5) 

See for example Appendix C in Ref. fzzll or Eqs. (6,7) 
in Ref. 12311 for more general formulae, allowing for non- 
equatorial measurements (ones with & z components). 

We have chosen the form of the Bell coefficients to coin- 
cide with the coordinates of an existing point in the set of 
quantum correlations because for this point this is the Bell 
inequality achieving the largest possible quantum value di- 
vided by the length of the vector defined by the coefficients. 
The reason we allowed only equatorial measurements to get 
the point in the correlation space was that maximum quan- 
tum value for the GHZ state can always be achieved for any 
full-tripartite correlation-type Bell inequality without tak- 
ing more general operators (see for instance 112311 ). 

We note that the quantum value we consider as above 
is usually not the maximum quantum value of the Bell in- 
equality. Even allowing only 3-qubit GHZ states there are 
usually quantum correlations giving larger quantum viola- 
tions. However, the improvement appears usually only in 
the fourth or fifth significant digit by the time the optimiza- 
tion procedure is finished, so it is unnecessary to determine 
the true maximum violation. The objective function has 



3 



TABLE I: Bell inequalities with critical visibility v or it < 0.5. The two cases with three measurement settings for two parties and one for 
one party, and the best ten cases with three measurement settings for two parties and one for the other two. Symbols C and Mijk denote 
the local bound and the Bell coefficients, respectively. 
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to be evaluated many times during the optimization pro- 
cedure, and just taking the sum of the squares of the Bell 
coefficients as the quantum value is much faster. In princi- 
ple, allowing other than 3-qubit GHZ states, possibly states 
in higher dimensional Hilbert spaces might further improve 
the quantum value, but we have found that for all tight Bell 
inequalities we arrived at in the second step this was not 
the case. In the second step, similarly as in the work of 
Ref. [25], we use linear programming to find the facet of 
the polytope in the direction of the vector defined by the 
Bell coefficients we got in the first step. In this step the 
critical visibility improves significantly most of the case. 
We note that usually we got the best final results not from 
the initial inequalities giving the lowest critical visibilities. 



Those most often led to inequalities that are equivalent to 
the Mermin inequality. 

IV. RESULTS 

For cases of no more than three measurement settings 
per party we have found no Bell inequality with critical 
visibility v cr u better than 0.5. For four measurement set- 
tings for one of the parties and three for the other two we 
have found two such inequivalent inequalities, V343 with 
VcrU = 0.49967, and V^ 3 with v crlt = 0.49972, see Ta- 
ble U We generated 30000 inequalities with the simplex 
downhill method starting from different random initial pa- 
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TABLE II: The best ten Bell inequalities with four measurement settings per party with critical visibility v cr u < 0.5. Symbols £ and 
Mijk denote the local bound and the Bell coefficients, respectively. 
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rameter values. In the second step we got inequalities equiv- 
alent with I/343 374 times, and with V343 68 times, respec- 
tively. The rest has v cr n = 0.5, many of them equivalent 
to the Mermin inequality, or v cr i t > 0.5. We have done the 
same number of attempts for three measurement settings 
for one of the parties and four for the other two. We arrived 
1372 times at inequalities with v cr u < 0.5, among them 
there were 70 independent ones, one of them equivalent 
with V343. Many of the inequalities were found only once, 
which indicates that there are certainly a lot we missed. The 
ten best ones we have got are shown in Table |U 



For four measurement settings per party we have made 
10000 attempts. In about one third of these attempts we got 
inequalities with v cr it < 0.5, of which over one thousand 
were inequivalent. The ten with the best critical visibili- 
ties are shown in Table HH It is interesting, that the integer- 
valued Bell coefficients are quite small for all of them. In 
Table HIl we show five more inequalities with four measure- 
ment settings per party with special properties. Two of them 
have only Bell coefficients and ±1, while the other three 
inequalities have coefficients symmetric in the exchange of 
all three parties. 
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TABLE III: Some special Bell inequalities with four measurement settings per party with critical visibility v cr u < 0.5. For two cases the 
absolute value of the Bell coefficients are no more than one, while three cases are symmetric for all participants. Symbols C and Myfc 
denote the local bound and the Bell coefficients, respectively. 
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34 
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-8 


34 


-15 
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-26 


-21 


19 


44 


-16 


-8 


-26 


-16 


-12 






-43 


-21 


-1 


35 


-21 


19 


44 


-16 


-1 


44 


-11 


-6 


35 


-16 


-6 
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-30 


-8 


35 


-7 


-8 


-26 


-16 


-12 


35 


-16 


-6 


17 


-7 


-12 


17 


44 



We have also made 10000 calculations for five measure- 
ment settings per party. In almost three quarters of the trials 
we got result with v cr a < 0.5. The best inequality, and 
some more inequalities have been found several times, but 
we have arrived at the majority of them only once. Some 
examples are shown in Tables [IV] and [V] The improve- 
ment in v cr it compared to the four measurement settings 
per party case is not very impressive. We show the three 
best cases, the second best in Table [Vl separately, for typo- 
graphic reason: the integer Bell coefficients are fairly large 
for that case. In TablellVlwe also show some special cases. 
V 5 55 is symmetric in all participants (the analytical solu- 
tion of this inequality will be given in Sec. [Vj, while 
and V555 have only coefficients and ±1. In a few cases 
we got inequalities equivalent to ones with less than five 
measurement settings per party. One example is V^ 55 , with 
only four measurement settings for Alice, for which v cr n is 
only very marginally worse than for V$ 55 . V454 has only 
and ±1 Bell coefficients and two participants have 4 mea- 
surement settings, only one has five. The total number of 
measurement settings is also 13 for V553, which also has 
Bell coefficients of absolute values at most one, and which 
is symmetric in the coefficients of Alice and Bob, who have 
five measurement settings each, while Cecil has only three. 
For this case v cr i t is quite small as well. 

For the critical visibility, we calculated the quantum 
bound with von Neumann projective measurements per- 
formed on components of a shared 3-qubit GHZ state with 
a see-saw method, similar to Ref. B24I1 . For all cases such 



calculation has given the exact quantum bound. This is true 
not only for the cases we reported here, but for all tight 
three-party full-correlation type Bell inequalities (nonzero 
coefficients only for three party correlators) we checked, in- 
cluding the ones with worse critical visibilities. We checked 
all cases we got with no more than four measurement set- 
tings per party, and also several larger cases, including ones 
with five measurement settings per party. We have done 
that by calculating the upper bound for the quantum value 
using semidefinite programming with the method invented 
by Navascues, Pironio and Acfn (NPA) l26ll . For the calcu- 
lation we used the code CSDP of Brian Bore hers $2~fo . In all 
cases we found the upper bound coincides with the result of 
direct calculation, which is actually a lower bound. In ap- 
plying the NPA method, it was enough to do it at level one 
+ab + ac + bc + abc in all cases. The notion of levels, and 
their notation is explained in Ref. l28ll . and also in Ref i29ll . 



V. STUDY OF A BELL INEQUALITY PROVIDING 
VISIBILITY LESS THAN 1/2. 

We consider a Bell scenario where each of the three par- 
ties can choose between five possible binary measurements. 
We denote the outcome of measurement j = 0, 1, . . . , 4 for 
Alice by Aj (Bj and Cj for Bob and Cecil respectively). 
Here we shall focus on the only symmetric inequality we 
found for this scenario (and can be found in Table IIVI ). 
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TABLE IV: Examples of Bell inequalities with critical visibility v C rit < 0.5 with five measurement settings for at least one party. Symbols 
L and My/t denote the local bound and the Bell coefficients, respectively. 
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which is given by the expression: 

VH =sym[-A B C 2 - A B C 3 + A B 2 C 4 

+ A Q B 3 d + A X B X C 2 - A X B X C Z - A 1 B 2 C 2 
+ A1B3C3 — 2A 2 B 2 C 2 + A 2 B 2 Cs — A 2 B^Cs 
+ A 2 B 4 C 4 + 2A 3 B 3 C 3 + A 3 B 4 C 4 }, (6) 

where the designation sym[X] indicates that the expres- 
sion X is symmetrized over all the parties, for instance 
sym[A B C 2 } = A B C 2 + A B 2 C + A 2 B C - In this 
5x5x5 Bell inequality the local bound is £ = 12 (all local 
probability distributions satisfy the inequality V^fg < 12). 
We now show that by performing local measurements on 



a 3-qubit GHZ state, 

\GHZ) = (|000) + |H1))/V2, (7) 

it is possible to obtain the value of Q — 24. 1699 for the Bell 
expression (0 above. This implies the visibility of v cr u = 
12/24.1699 = 0.496485 showed in TablellVl 

The local measurements used by each party are of the 
simple form Aj = cos (pj<r x + sin ipj& v (with j = 
0, 1, . . . , 4), where a x and a y are Pauli matrices. Moreover, 
the five local measurements are the same for each party, i.e. 
Aj = Bj = Cj, j = 0, 1, 2, 3, 4. In this case the 3-party 
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TABLE V: The Bell coefficients of V555, the Bell inequality with five measurement settings per party with the critical visibility v cr it = 
0.496059, the second best one we found. The local bound is C = 13180. 

-271-124 -418 522 251 922 355 8 522 -37 540 206 69 90-641 -294 437 853 90-212 65 710 210 -565 

-514 251 736 596-403 -138 -45 -249-612 736 -217-641-350 932 853-212 -45 82-514 418-557-400 -98-1277 

-446 271-1146-457 514 -636-922 -482 612 130 546-540-996-167 209 -347 294-477 175-853 775 -65-145 147 -418 

-457-522 661 0-596 604-514 858-336 612 -159 -90-495 426 159 -90 763 604 -90 147 -8-465-514 90 

1146-418 669-661 736 466 8-1915-858-233 996 69 368 503-342 461 853 342-771 -37 145 210-510 465 -400 



correlation term has the following form (see Eq. ©), 

AiBjCk = cos((fi + cpj + ip k ). (8) 

The optimal measurements are given by cpo = 5tt/8, (fx = 
7r/2, ipi = 7r/8, and ip2 + </?3 = 37r. Leaving ip2 as a free 
parameter, we obtain the following closed formula for the 
Bell value, 

S = — 12 cos 3<^2 + 12 cos if2 (sin ip2 — 1) — 12\/2sm<^2. 

(9) 

Our task is now to perform optimization over the angle (f2- 
By inserting tp 2 = 3.73842 into (0, we obtain the value of 
Q = 24.1699 for the expression ©■ 



VI. CONCLUSION 

The non-local properties of the noisy three-qubit 
Greenberger-Horne-Zelinger (GHZ) states (defined by 
Eq. parameterized by the visibility < v < 1 were 
investigated. Based on the violation of the 2 x 2 x 2-setting 
Mermin inequality, p v is non-local for the parameter range 
1/2 < v < 1. In this study we presented several Bell 
inequalities which beat the limit of v = 1/2, thereby an- 
swering the question raised by Kaszlikowski et al. lfl3ll . 
In particular, the lowest threshold visibility we found is 
v = 0.496057, attainable with 5x5x5 settings, whereas 
the most economical one overcoming v = 1/2 corresponds 
to 3 x 3 x 4 settings. The method which enabled us to ob- 
tain these results, and in particular the about 10000 tight 
Bell inequalities going below v — 1/2 were also discussed 
in detail. 

We expect that the presented method in Sec.|III]would be 
efficient to test the non-locality of other states as well. For 



instance, in the two-party 2-qubit scenario it might prove 
useful to find Werner states with better visibility thresholds 
than the presently known ones iTToh (or to find ones which 
overcome the 1/V2 limit corresponding to the CHSH in- 
equality with fewer number of settings). Our technique may 
also be suitable to generate Bell inequalities for 3-qubit W 
states or 4-qubit Dicke states Bill , cluster states JH] 
with improved noise resistance. 

Further, it would be interesting to study whether the tech- 
nique presented could be applied to find Bell violation with 
lower detection efficiencies than the presently known ones 
both for bipartite and multipartite settings. In the bipartite 
case, according to Ref. 0311 . one may postulate the d x d 
state to be = e|00) + |11) + . . . + \d - 1, d - 1) (here 
the state is written in an unnormalized form), e is a small pa- 
rameter and the measurements are rank-1 projectors. Then, 
by using appropriate parametrization, one might carry out 
the heuristic search of Sec. [Ill] for 2-outcome Bell inequal- 
ities to lower the detection efficiency required for closing 
the detection loophole. 

In the tripartite case, it is known that for the GHZ state 
the minimum detection efficiency based on the Mermin in- 
equality is 75% ll34ll . On the other hand, the best protocol 
which simulates GHZ correlations (arising from equatorial 
von Neumann measurements) with detection efficiencies of 
50% is due to the recent work of Ref. 1 1 611 . It would be inter- 
esting to decrease thegap by considering inequalities with 
more than 2 settings 113511 in conjunction with the technique 
used in this paper. 
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